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The origin of large scale cosmological magnetic fields remains a mystery, despite the continuous 
efforts devoted to that problem. We present a new model of magnetic field generation, based on 
local charge separation provided by an anisotropic and inhomogeneous radiation pressure. In the 
cosmological context, the processes we explore take place at the epoch of the reionisation of the 
Universe. Under simple assumptions, we obtain results (i) in terms of the order of magnitude of 
the field generated at large scales and (ii) in terms of its power spectrum. The amplitudes obtained 
(B ~ 8. 10~ 6 uG) are considerably higher than those obtained in usual magnetogenesis models 
and provide suitable seeds for amplification by adiabatic collapse and/or dynamo during structure 
formation. 

PACS numbers: 98.80.-k; 98.65.-r; 95.30.Qd 



I. INTRODUCTION 

According to observations, magnetic fields are pervad- 
ing the whole Universe. In particular, both synchrotron 
emission and Faraday rotation measurements reveal a 
magnetic field present on large scales. On galactic scales 
[1|, on galaxy cluster scales S and even on larger scales 
[3|, the observed magnetic field strengths range from 0.1 
micro-Gauss to a few tens of micro-Gauss. 

The origin of the magnetic field is, however, still a mys- 
tery. Clearly, the fields cannot be generated directly with 
the observed amplitudes since, in that case, they would 
strongly affect the formation of structures jjj. Therefore, 
the fields must arise in the form of weak seeds which are 
subsequently, in some way, strongly amplified after the 
onset of structure formation. Various mechanisms have 
been proposed for the generation of such seeds. Basically, 
they can be sorted in two types. 

In cosmology, the mechanisms of the first type oper- 
ate before matter-radiation decoupling. They take ad- 
vantage of high energy physics processes such as infla- 
tion JE], ^, or phase transitions in the early Universe 
H ||. However, these mechanisms are still unsatisfactory 
in two respects. First, they do not provide a real predic- 
tion since the generated fields can be as weak as 1CU 65 
Gauss and as strong as 10~ 9 Gauss. Moreover, a prop- 
erty shared by fields created before decoupling is that the 
stronger their amplitude, the smaller the scale on which 
they are generated (see Battaner & Lesch |pxj| , and ref- 
erences therein). Furthermore, under the assumption of 
frozen magnetic flux, these fields decay following cosmo- 
logical expansion as (1 + z) 2 . Therefore, fields of primor- 
dial origins (generated before nucleosynthesis, say) are 
extremely weak on galactic scales and are very strong on 
much smaller scales. Finally, extremely tight constraints 
on pre-BBN magnetic fields from gravity wave produc- 
tion have been recently derived pi, possibly ruling out 



a majority of the proposed models. 

The mechanisms of the second type produce magnetic 
fields after matter-radiation decoupling. They rely essen- 
tially on the battery effect jl^] which was first introduced 
to explain stellar magnetic fields and is basically a conse- 
quence of charge separation. Several authors explored 
the same mechanism in various cosmological contexts 
]l3| , |l5| . The intensity of the fields generated is found 
to be roughly of the order of 10 -20 Gauss. These weak 
seeds are then amplified on galactic scales by a dynamo 
sustained by turbulence and differential rotation in pro- 
togalaxies. As noticed by Grasso & Rubinstein Q, this 
type of mechanism can hardly be responsible for mag- 
netic fields on galaxy cluster scales. Consequently, we 
would have to deal with different mechanisms to explain 
the origin of fields in galaxies and clusters, which exhibit 
similar properties. Moreover, this mechanism cannot re- 
ally account for the strong fields observed in objects at 
high redshifts where galaxies did not have enough time 
to rotate and amplify the seeds efficiently. 

Finally, whatever the problems of each type of mecha- 
nisms, the seeds produced are generally small. As such, 
they require the action of a strong process capable of 
amplifying the magnetic amplitudes up to the values ob- 
served in galaxies. The dynamo mechanism evoked above 
has long been considered suitable for that purpose but it 
is, however, subject to controversy nowadays. In fact, 
it has recently been pointed out jl7| that dynamo tends 
to amplify the small scale magnetic fields in a more effi- 
cient way than fields on large scales. The consequence of 
this is that small scale fields reach rapidly dynamically 
important strengths and stop, in back-reaction, the dy- 
namo mechanism itself. This happens long before large 
scale fields can reach the observed equipartition values. 

Significant progress has been achieved very recently 
and the back-reaction of small scale fields is now bet- 
ter understood. In particular, the importance of mag- 
netic helicity conservation in dynamo quenching has been 
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demonstrated in ideal cases [18 



However, for 



their application to real astrophysical systems such as 
galaxies, these theories require further developments. 
Small scale fields thus remain a possible threat to dy- 
namo amplification mechanisms. 

In this paper, we propose a mechanism which appears 
to be relatively free of the problems outlined above. Ap- 
plied to cosmology, it operates after decoupling in the 
natural context of inhomogeneous reionisation, for a red- 
shift z between 6 and 7 [|lj [2^] . Similarly to the battery 
effect, it relies upon local charge separation provided in 
our case by radiation pressure. The physics involved is 
fairly simple and was nevertheless left aside until now, 
surprisingly enough. Notice however that Subramanian 
et. al. ]23| proposed a mechanism (which has been ex- 
plored numerically by Gnedin et. al. |Q) for the gen- 
eration of magnetic field seeds at reionisation, and ob- 
tained magnetic fields of the order of 10 -20 — 10~ 19 Gauss 
on galactic scales. Their solution for magnetogenesis re- 
lies explicitly on the actual battery effect which requires 
non-parallel temperature and electron density gradients. 
While exploring the physics at a comparable epoch, we 
consider here a different and more efficient process, as 
will appear below. 

In the following section, we describe in details the 
mechanism of our model. Then, in Sect. 3, we present 
the results obtained for the cosmological case, emphasis- 
ing on the magnetic field power spectrum and its order of 
magnitude at protogalactic and galactic scales. We also 
apply our mechanism to the case of interstellar magnetic 
fields. Finally, we discuss the model and its limitations 
in Sect. 4. 



II. PHYSICAL PROCESSES 
A. Description 

Consider a radiation flux propagating through an 
ionised plasma. Far away enough from the dominant 
source, the radiation is anisotropic and its propagation 
direction defines a local preferred axis. Since the scatter- 
ing cross sections of electrons and ions are very different 
(ctt cx mass -2 ), the radiation pressure combined to the 
anisotropy of the flux lead to a charge separation, which 
in turn creates uniform electric fields on large scales com- 
pensating the difference in radiation pressure. 

Furthermore, the radiation propagates through a 
medium where the matter is distributed inhomoge- 
neously and some of the over-densities act as screens 
for the radiation. Moreover, the sources of the radia- 
tion are not distributed homogeneously. Thus, because 
of the source distribution and the gas density inhomo- 
geneities, the radiation flux and particularly its ionising 
part, is anisotropic and inhomogeneous far away from 
the sources. To first order, the inhomogeneities of the 
flux remain unchanged along the propagation axis over 
regions which correspond to typical distances between 
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FIG. 1: Schematic view of the magnetic field generation pro- 
cess. The source emits a radiation flux $ which, at some 
distance, defines a preferred direction (Oz) (anisotropy). It 
ionises the diffuse gas, and the radiation pressure locally sep- 
arates electrons from ions. A static uniform electric field E 
is generated in response to <1? and the electrostatic force bal- 
ances the radiation pressure. The situation is different behind 
an overdensity: In the shadowed region, the radiation flux is 
weaker (<l?i <C <&), there is no static equilibrium and free elec- 
trons are set into motion (v B ). The electric current j created 
this way closes up outside of the shadowed tube and induces 
a magnetic field B on the scale of the overdensity. 



matter density fluctuations. In such regions, the vari- 
ations in the flux allow electrons to move and electric 
currents to develop. The latter finally induce a magnetic 
field on the length scales over which the radiation flux 
varies (see figure [j] for an illustrative description of the 
processes as they would work for one single source and 
one single density inhomogeneity) . 

The situation described above is likely to exist at the 
time of the reionisation of the Universe when the very 
first sources (population III stars or first galaxies and 
quasars) turn on and ionise the neutral primordial gas, 
preferentially in the low density regions. Collapsing or 
collapsed structures play the role of the matter inhomo- 
geneities leaving their imprints in the radiation flux. 



B. Formalism 

We study the situation in a region far enough from the 
dominant ionising source where the direction of the ion- 
ising flux «I> can be considered constant. This direction 
defines the axis (Oz) (sec fig. [l]). The anisotropy of the 
flux is essential in that it leads to a large scale charge 
separation which in turn generates a large scale electric 
field. We could have equivalently considered a globally 
isotropic flux with an additional anisotropic component. 

Further, in our model, we describe the anisotropic ion- 
ising flux $ as a homogeneous background <& with an 
inhomogeneous part superimposed on it, 



3> = 3>o + 3?i 



(1) 
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where 3?i accounts for the inhomogeneities which are 
taken to be small with respect to the constant part 
4>o = (*}, namely 3>i = f(r) <fr with f(r) <C 1. As 
we wrote above, the region which we study is far enough 
from the dominant source. At such a distance, the den- 
sity fluctuations in the medium along the line of prop- 
agation (Oz) of the radiation have left imprints on the 
ionising flux which is therefore inhomogeneous. In other 
words, the 3>i part of the flux takes into account the 
inhomogeneities due to both the distribution of sources 
and the variations of the flux along the line of propa- 
gation (Oz) due to the variations in the matter density. 
In the regions of study, the properties of the medium 
are supposed not to vary significantly along the line of 
propagation. Provided this assumption, the spatial de- 
pendency of $1 is confined to the plane (xOy), which 
implies that derivatives with respect to z can be taken 
equal to zero, and f(r) — f(x,y). Finally, we also con- 
sider the stationary regime, and therefore all the time 
derivatives are nil as well. 



C. Maxwell equations 

Together with the assumptions specified above, the 
usual Maxwell equations 
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(p and j are the charge and the current densities respec- 
tively, c being the speed of light) reduce to the following 
simple relations: 
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for the Faraday equation (2b) 
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for the Ampere equation (gd|) . Relations (g) indicate that 
the electric currents constitute the only sources of mag- 
netic field in the problem. We consider that the medium 
is globally neutral. The charge balance between free elec- 
trons, of number density n c and charge q e , and ions, of 



density n\ and charge cfi, implies n\q\ + n c q c = 0. The 
electric current is j = n c q c v c + niqiVi = n c q c (v c — v{). In 
the following, we suppose, as customary in similar prob- 
lems, that the ion velocity is negligible with respect to 
the electron velocity, and we write therefore j — n c q c v . 

Combining the Coulomb and the Ampere relations 
with the dt = 0, d z = conditions, the charge conserva- 
tion equation takes the following form: 



d x j x + d y j y = 0. 



(7) 



Within the stationary assumption, if we take the curl of 
equation (Eq) we obtain 



V 2 B = -^Vxj 



(8) 



which, because of our assumptions on the spatial depen- 
dency of 3>i, reduces to the following relations between 
the magnetic field and the current density: 



v 2 s, =-— a,, 
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47T 



V 2 By = —3 xJz 

c 

47T 

V 2 B Z = —{dyj X - d X jy) 



(9a) 
(9b) 
(9c) 



We will later combine the latter relations with the equa- 
tion of motion. 



D. Equation of motion 

Given the relation outlined in the previous section be- 
tween the current j and the electron velocity v c , the 
equation of motion for the electrons can be re-written 
in terms of the current evolution: 



— -£ = q c E+ J x B —j + F cxt 10 

n e q e dt n c c n c q c 

where m e is the electron mass, v c is the mean collision 
frequency associated with momentum transfer between 
electrons and ions, and _F cxt stands for the external forces. 
The derivative on the left hand side of this equation is the 
whole Lagrangian derivative, i.e. 4z = d t + v ■ V which, 
in the stationary regime, reduces to ^ — v ■ V. On the 
right hand side of equation (]l0|), the first two terms are 
the Lorentz force, while the third term represents the 
Ohmic current dissipation which depends on the physi- 
cal state of the medium. For a partially ionised gas, we 
should add to ( ^p| ) a term representing the force exerted 
by the neutrals on the electrons. However, we concen- 
trate on the case where the medium is ionised enough for 
the collisions with the neutrals to be negligible. Thus, we 
do not consider here this effect. Neglecting also gravity, 
the only external force applied on the electrons is the ra- 
diation pressure force, which in terms of the ionising flux 
reads F cxt = — &t * {hv is the mean energy of ionising 
photons and ctt is the Thomson scattering cross section) . 
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E. Linearised problem 

Assuming the inhomogcncitics f(x,y) of the ionising 
flux to be small, we can linearise equation (JT"o|) - Identi- 
fying the different terms of the same order leads to the 



following relations: 




m c v c 

*^X 9 Jx 

n c q£ 


(11a) 


_ m c v c . 

V 9 Jv 
n c q 2 c 


(lib) 


i hv 
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(11c) 


and 





. n c q c hv 

JA x ,y) = °T $o j(x,y) 

m c i> c c 



(12) 



where second order terms have been neglected. 

Equation (11c) shows simply the equilibrium between 



the radiation pressure and the electrostatic force exerted 
on the free electrons, whereas equation (|l2| ) relates the 
induced current to the flux inhomogeneities. The current 
changes direction depending on the sign of the deviations 
fix, y) from the mean value of the flux, and therefore the 
current lines are closed over scales similar to those of the 
flux variations in the (xOy) plane and comparable to the 
coherent scale of the anisotropy along z. 

Taking the curl of equation ( |l0| ) and keeping again only 
the first order terms gives 



c c V m c v c 



(13) 



where equation (^|) has been used to eliminate the cur- 
rent. Again, due to the particular g eome try of our prob- 
lem and taking profit of equations (11a) and ( lib ), the 
equation above can be rewritten as 



V 2 B X = -d v [Kf(x,y)} 

«/0,y)] 



^ By 
V 2 B Z 



= d x 
= 



(14a) 
(14b) 
(14c) 



with k = *2L J2^k- hiL cr T $ . The generated magnetic 
field thus depends on the physical state of the medium 
through the friction coefficient . It depends also on 
the flux through both its mean intensity $o and its spa- 
tial fluctuations f(x,y). 



III. RESULTS 



The solutions to (14) give the generated magnetic field. 
Formally, we obtain 



B x (r) 
B v (r) 



hi\r'-r\ d y («/(r)) dV 



ln|r'-r| d x (nf(r)) dV 



(15a) 
(15b) 



Of course, it is also poss ible to find a solution for B z . In- 
tegrating equation (14c) would require precise boundary 
conditions. But in the model presented here, there are a 
priori no such obvious physical conditions, and imposing 
for instance B z = at some distance would introduce 
an arbitrary non physical scale. Actually, recall that the 
study we present here applies only to the domain where 
an axis (Oz) is defined by the anisotropic flux and in 
which all the z-dependencies can be projected onto the 
(xOy) plane. In reality, many zones with locally con- 
stant $ may coexist and the local directions of the flux 
would be independent of each other. At the common 
borders of these regions, the three components of the 
field would get mixed. In that case, the only pertinent 
boundary conditions for B z would be fixed by the x and 
y components of the fields generated in similar surround- 
ing regions. Given the restrictions provided by our initial 
assumptions, including such a complex description of the 
situation is obviously beyond the scope of this paper. 

Therefore, our intention here is not to find the exact 
value of the magnetic field at each point in space but 
rather to address the question of the origin of large scale 
magnetic fields. For that purpose, provided the relations 
found in ( |l4| ) , it is more interesting to derive global prop- 
erties of the magnetic field such as its order of magnitude 
and its power spectrum. 



A. Order of magnitude 

In the previous section, we found the equations which 
relate the transverse components of the magnetic field to 
the ionising flux perturbations (Eqs. |lj) . From those 
relations, we can estimate the order of magnitude of the 
generated field. Relations ([L41) can be expressed as 



B^HifR 



(16) 



where R is the characteristic scale over which / (hence 
B) varies. 

As we wrote earlier, the friction coefficient = 
appearing in k, depends on the physical state of the 
plasma. In particular, it varies according to the am- 
plitude of the magnetic field, and two regimes can be 
distinguished. When the electron-ion collision frequency 
v cl is greater than the cyclotron frequency cj c , we are in 
the non-saturated regime, v c = v B1 and the electrical con- 
ductivity is given by a c = er w 6.5 10 6 T 3 / 2 e.s.u [g5| for 
a fully ionised medium. But as soon as the magnetic field 
is strong enough to have lo c 3> v c \ (i.e. in the saturated 
regime), the actual electron-ion collision frequency, thus 
also the conductivity, are modified anisotropically. When 
the particle radius of gyration becomes smaller than the 
mean free path, the situation is unaffected along the mag- 
netic field lines whereas collisions transverse to the lines 
are less likely to occur. Therefore, to the first order, the 
conductivity along the magnetic field lines is unchanged, 
whereas in the transverse direction (along the axis (Oz) 
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in our case), a c takes the form a c ~ <tq i (known as 

' ^ci 

the Pederson conductivity). 

Finally, inserting these expressions into ( |l6| ) and 
rewriting hv^o as j^jyi (L is the luminosity of the ionis- 
ing source and D the distance to it), we end up with the 
following order of magnitude estimates: 



B 



3.09 10" 
R 



f 



T 



10 4 K 

LD- 2 



3/2 



100 kpc 10- 8 W.m- 



Gauss 



(17) 



for the magnetic field seed generated in a fully ionised 
plasma in the non-saturated regime, and 



B 



3.14 10~ 2 / 



R 



1/3 



( T 



VlO 4 K 
LD- 2 



1/2 



100 kpc 10" 8 W.m" 2 



1/3 



5 S 2/3 Gauss (18) 



in the saturated regime (B s = 2Ttv e - 1 m e ce~ 1 is the satu- 
ration magnetic field which depends, through u e i, on the 
density and the temperature of the medium) . 



B. Application to cosmological magnetic fields 

We now want to apply the result obtained above to the 
case of cosmological magnetic fields. Since the plasma 
which we consider is (almost) fully ionised at z ~ 7, we 
take the temperature T to be of the order of 10 4 K. For 
a flat cosmology in a matter dominated Universe with 
fl b hl = 0.020 §§, the gas density is n(z) = 3.8 10" 7 (1 + 



3 and the corresponding saturation magnetic field 
■ 1.4 lO" 11 T~ 3 / 2 (1 + zf Gauss which, for a 



K, gives B s ~ 1.4 10~ 17 (1 + z) 



z) cm" 
is B s r 

temperature of 10 4 
Gauss. 

The term LD~ 2 can be estimated by calculating the 
number of photons necessary to reionise the Universe. 
Basically, at the reionisation rcdshift z- loni the integrated 
photon density production must be equal to the gas den- 
sity: 



CN s (z) dt = n(zum) 



(19) 



where C is the number density of photons emitted per 
unit time and N s (z) is the density of emitting sources 
per logarithmic interval of luminosity at redshift z. The 
major contribution to reionisation comes from photons 
emitted roughly at redshift z- lon , so that the above rela- 
tion can be rewritten as 



^-^s^iorO^u^ion) — ^-(^ion)- 



(20) 



In a flat, matter dominated cosmology, the age of the 
Universe is t u (z) = jHq 1 (1 + z)~ 3 / 2 . Therefore, 



CN s (z ion ) = 5.7 10- 7 ff (l + ^ion) 



9/2 



(21) 



The term LD 2 is then written 
LD~ 2 = 8.1 10~ n /in /3 



20 eV 



L 



1/3 



(1 



2/3 



3 W.m" 



(22) 



,_t _i is the reduced Hubble con- 



where hn = -nr n , , . 

u 100 kpc.Mpc 

stant and hv is the photon energy (in units of 20 eV). 
Taking a luminous quasar as the nearest ionising source 
(responsible for the anisotropy of the flux), L ss 1O 12 L0. 
With hv w 20 eV, and considering ho « 0.70, equation 
(pi) gives 



LD' 



i 6.38 10- 3 (1 + z ioa f 



(23) 



10~ 8 W.m- 2 

At a redshift z- lon ~ 7, this corresponds to D m 3.12 
Mpc. Notice that this is an underestimate of the mean 
distance to the nearest quasar since we assumed implic- 
itly in equation (|l^) that luminous quasars were alone 
responsible for the ionisation of the medium. 

With the values of the density and the temperature 
reached at high redshifts, it is easy to check that the es- 
timated non-saturated magnetic field is greater than the 
saturation value B s . Therefore, to obtain the correct or- 
der of magnitude for B, we have to consider the saturated 
regime. Inserting the figures derived above into equation 
(|l8|), we obtain 



B - 1.56 10" 



R 



100 kpc 



1/3 



(1 + z ion ) 3 Gauss (24) 



where we have taken T w 10 4 K and the degree of fluc- 
tuations in the ionising flux / ~ 10%. 

The length R represents the typical scale of the inho- 
mogeneities in the radiation flux. As we already men- 
tioned, the flux inhomogeneities take into account simul- 
taneously both the density fluctuations in the medium 
and the variations due to the spatial distribution of ion- 
ising sources. Therefore, the scale R in our calculations 
is not easy to define precisely without proper statisti- 
cal estimates of the source and gas density distributions. 
Nevertheless, we can at least estimate reasonable bounds 
for it. On large scales, the mechanism proposed here 
remains valid as long as the flux holds its property of 
anisotropy. Clearly, on scales larger than the mean dis- 
tance D between the considered ionising sources, this is 
no more the case, and therefore D defines approximately 
an upper limit for R. On small scales, the variations of 
the flux would be due essentially to the matter density 
fluctuations. More precisely, the lower limit for R would 
be set by the typical size of the smallest halos which 
could survive photoevaporation by the ionising flux it- 
self. However, for the simple order of magnitude eval- 
uation intended here, it is sufficient to choose R as the 
size of a protogalaxy, R w 100 kpc, in which case, the 
magnetic field seed at protogalactic scales is 



B ~ 8. 10~ 12 Gauss. 



(25) 
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This result is very interesting since it is several orders 
of magnitude higher than fields found in usual plasma 
physics magnetogenesis models. 

One may ask whether turbulence will affect the mag- 
netic field generation mechanism. Indeed, in the vicinity 
of the ionising source, the gas motion is highly turbulent, 
but in fact, the actual effect of turbulence is not easily 
predicted. Presumably, turbulent motions will strongly 
fold the field lines on small scales leading to local mag- 
netic field amplification and/or to magnetic reconnection. 
Moreover, turbulence may even prevent electric currents 
from settling over large distances, thus inhibiting the 
mechanism presented here. 

However, even in presence of turbulence, the ionising 
radiation will be able to escape the vicinity of the source. 
Thus, far enough from the source, which is the case of in- 
terest explored here, electric currents will be able develop 
and induce magnetic fields. On large scales, the growth 
of the density fluctuation field is in the linear regime and 
the collapse of a large scale structure is still in its early 
stages where the fluid motion is not turbulent and rel- 
ative displacements are small. Nevertheless, given the 
relatively large amplitude derived in Eq. [24|, we might 
want to check that the magnetic pressure will not affect 
the collapse of the structure. The velocity 5vk of gas 
within a density perturbation 5k on a scale R = 2ir/k is 
given by 



\5v k \ ~ ?-8 k ~ 11.1 4(1 + *) 3/2 tt?- km.s- 1 
K 1 Mpc 



(26) 



in the linear regime and for a flat, matter dominated 
Universe with Hq ~ 70km.s _1 .Mpc~ 1 . Thus, the kinetic 
energy density is written 



1 2 
E K = -Pb<>k\5vk\ 

E K w 2.3 lCr 19 (5j* 



(27a) 



R 



k \ 1 Mpc J 
x(l + z) 6 erg. cm -3 , 



where pb ~ 3.76 10 31 (1 + z) 3 g.cm 3 is the mean baryon 
density, whereas the magnetic energy density is 



Eb = 

8tt 



(28a) 



R 



2/3 



e b « 2.1 icr 28 / 2 / 3 , 

1 1 Mpc J 

x(l + z) 6 erg.cm~ 3 , (28b) 
from Eq. ^J. The ratio of the two energy densities, 

-4/3 



Eb 
Ek~ 



9.1X0- 10 / 2 /3*-3 ( -J*—) 
k V 1 M P C / 



(29) 



does not depend explicitely on redshift. At a scale 
R ~ IMpc and considering / ~ 5 k ~ 0.1, this gives 
Eb/Ek ~ 2. 10 -7 . According to the last relation above, 
Eb and Ek become comparable on scales smaller than a 



parsec. At such small scales, structures are deep in the 
non-linear regime which we do not consider here. There- 
fore, the pressure provided by the magnetic field gener- 
ated by our mechanism on large scales is well below the 
kinetic energy density. This guarantees that the collapse 
of the structure will not be modified by the presence of 
the magnetic field. Thus, as stated above, as long as 
the structure is in the linear regime, the collapse will not 
affect our mechanism. 

On the contrary, the collapsing gas will drag the field 
lines and amplify the field strength eventually. In fact, 
the amplification by adiabatic collapse gives a galactic 
magnetic field sufficiently high to serve as seed for dy- 
namo. The galactic value obtained from (^5|) after adi- 
abatic collapse (and before further contraction by radia- 
tive cooling) is 



R_, = 5 2/3 B ~ 2.7 10" 



Gauss 



(30) 



(where 5 C ~ 200) from which dynamo amplification needs 
less time than in usual models in order to reach the fields 
of a few micro-Gauss observed in galaxies. 



C. Magnetic field power spectrum 

Let us start from the magnetic field two-dimensional 
spatial correlation function 



£(r) = (B(r)-B{r + r')) 



(31) 



where r and r' are confined to the (xOy) plane. Inserting 
the Fourier transform of B in that expression, we end up 
with 

£(r) = (^) f B fe B_ fe e- ife M 2 fe 



(27b) f(r) 



\B k \ 2 e- lkr d 2 k 



(32) 



where the last line comes from the reality of B. The 
quantity we are now interested in is the power spectrum 
of the magnetic field Pb(&), i.e. the angular average of 
the spectral correlation function £(fc): 



(33) 



2tt Jo 



where £(fe) is the Fourier transform of £(r). It is straight- 
forward to see that £(fc) = \Bk\ so that the power spec- 
trum is written 



^ Jo 



(34) 



Restricting ourselves to the case of the non-saturated 
regime, we can consider the conductivity to be constant 
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and we can replace a e b y Co i n n. If we then Fourier 
transform equations (14a) and ( |l4b| ), we get 



^ 2 — k k y fk 
k = k k x fk 

which finally allows us to obtain 



;\fk 



k 2 



(35a) 
(35b) 



(36) 



Thus, the mechanism we propose generates magnetic 
fields preferably on large scales whereas small scale fields 
are partially suppressed. This is very interesting with 
respect to dynamo theories. Since in our model the small 
scale fields are partially damped, the problem of their 
amplification evoked in the introduction is less crucial. 
However, the exact expression of the magnetic field power 
spectrum Pe(fc) depends on the shape of the flux power 
spectrum \f k \ 2 . 

The fluctuations f(x,y) of the ionising flux are directly 
related to the fluctuations of the optical depth according 
to 



1 



where T\(x,y) — —ainh / 5(r)dz 
Jo 



(37) 
(38) 



(<7i is the ionisation cross-section, is the mean particle 
density and S represents the density fluctuations) and the 
fluctuation power spectrum can be expressed in terms of 
the power spectrum of the optical depth \r k \ ■ 

Following Bartelmann & Schneider [^7j in their deriva- 
tion of the Limber's formula in Fourier space, we relate 
the latter to the density power spectrum |(5fc| 2 as 



\Tk\ 



d 2 k' A 

2tt 



2 KV 
— — sin — — 

k'^l 2 



(39) 



The last integral in the equation above represents in 
Fourier space the projection on the (xOy) plane of the 
matter density over a distance I. Density fluctuations on 
longitudinal scales larger than I do not contribute to the 
projected power spectrum and the integration therefore 
must be carried over the interval ku € [^f-,+oo[. 

In the optically thin limit, Pb(^) oc |ffc| 2 /k 2 . Assum- 
ing a power law density spectrum \6k\ oc fc™, we calcu- 
lated the magnetic field power spectrum numerically in 
the non-saturated case. Results shown on figure || for the 
cases n = — 1 and n = — 2 assume the integration depth 
/ = 1 Mpc. The contribution of density fluctuations on 
scales larger than / to the optical depth fluctuations in 
the plane (xOy) is independent of the scale. Therefore 
\fk\ 2 is scale independent and -Pb(^) k~ 2 for k < 1/7. 
Note however that this part of the power spectrum is 
of little pertinence since it corresponds to scales larger 



10 



10 - 
10 



10 



-3 




FIG. 2: Power spectrum of the generated magnetic field in the 
non-saturated regime. A power law \St\ 2 oc k n power spec- 
trum has been assumed for the density fluctuations. Solid 
line: n — — 1; dashed line: n = —2. Normalisation is arbi- 
trary. Slopes are indicated by dotted lines. 



than I ~ 1 Mpc. These scales are comparable to, or 
larger than the mean distance between ionising sources 
and at those scales, our assumptions break down (lost 
flux anisotropy). On scales smaller than I, the magnetic 
field power spectrum is strongly decreasing. Depending 
on the spectral index n chosen for the density power spec- 
trum, small scale magnetic fields are strongly damped. 
For n w —1, corresponding roughly to the spectral in- 
dex on galaxy cluster scales, the magnetic field power 
spectrum goes like k~ 3 , whereas on galactic scales where 
n ~ —2, small scale magnetic fields are damped even 
more strongly, Ps(k) oc k~ A . This result makes our mag- 
netogenesis model a very interesting ingredient for dy- 
namo amplification theories for which small scale fields 
seem to be problematic, as explained in the introduction. 

We have nevertheless to bear in mind that the results 
above are derived in the non-saturated regime. A simi- 
lar study is currently undertaken for the saturated case. 
However, we might expect that even in the saturated 
case, small scale fields will be largely damped. 



D. Application to the Interstellar Medium 

Magnetic fields in galaxies are organised on large scale 
and have intensities of a few microgauss (see, e.g., Beck 
p8[ , and references therein). Observations show that the 
structure of the magnetic field in the Interstellar Medium 
(ISM) is correlated with the matter distribution comfort- 
ing the growing evidence that the ISM is structured not 
only by gravity but also by supersonic alfvenic turbu- 
lence. 

The situation presented in Sect. II A occurs naturally 
in the ISM, especially for clouds in the vicinity of OB 
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star associations. It is thus legitimate to ask whether 
our mechanism can contribute significantly to the ob- 
served magnetic field strengths and possibly account for 
the poloidal geometry of magnetic fields observed around 
molecular filaments [^9) , as suggested by figure [l]. 

As long as the electron-ion collisions dominate over 
electron- neutral collisions, we can apply directly the re- 
sults derived in Sect. Ill A. The electron-neutral collision 
cross section is of the order of a few 10~ 15 cm 2 |3(| . In the 
neutral diffuse ISM where T ~ 1. 10 4 K, n H ~ 1. cm" 3 
and — ~ 1. 10~ 2 , the electron- neutral and the electron- 
ion collision frequencies are v cn ~ 1.5 10~ 9 s" 1 and 
u e i w 9.6 10 -7 s _1 respectively. Therefore, the ionised 
component dominates in the diffuse neutral ISM. In in- 
terstellar clouds, where hydrogen is neutral but carbon 
is still ionised, T = 30 to 100 K, n B = 100 to 1000 cm" 3 



and ^ 
while v c 



~> 1. 10 , we have v t 
= 7. 10~ 4 to 3. 10 



-8 



cn 
-2 



= 1. 10~ 8 to 9. 10~" S" 

3 _1 and the contribution 



of ions is even more dominant. Thus, we can ignore the 
presence of neutrals and apply equations ( o) and (|lj) . 

When the electron-ion collisions dominate, the unsat- 
urated magnetic field is independent of the density and 
depends only on the temperature. Taking the photon 
flux given by a 10 6 L Q OB association at 10 pc and a 
transverse gradient created by a 1 pc interstellar cloud, 
we get 



B unsat = 1.25 10- 10 / T 3 / 2 Gauss 



(40) 



or -Bunsat ~ 1-25 lO -4 / Gauss in the diffuse medium 
and B unsa t ~ 1.5 lO -8 / Gauss in clouds. The critical 
fields -Bcrit where the Pederson conductivity saturates are 
2. 10- 10 to 1. 10- 8 and 3.5 10" 13 Gauss in the clouds and 
the diffuse medium respectively. We thus see that the 
saturation plays little role in clouds but reduces by a large 
factor the field in the diffuse medium. The saturated 
magnetic fields which are produced by our mechanism 



B. 



(-B unsat B 2 rit ) 1/3 



(41) 



that is -B S at = 2.4 10 10 J 1 / 3 Gauss in the diffuse medium 
and 2. 10~ 9 Z 1 / 3 to 1.5 10- 8 J 1 / 3 Gauss in clouds. For 
the specific case of molecular clouds, the induced mag- 
netic field is roughly one order of magnitude smaller than 
in clouds. We therefore conclude that in the context of 
the ISM, our mechanism is a subdominant source of mag- 
netic fields as compared with galactic dynamo and tur- 
bulence. In fact, the regime in the ISM is highly non 
linear, and the amplitude of relative motions of parti- 
cles is large. Therefore, the processes discussed in Sect. 
[II B dominate, and turbulence will possibly not allow 
electric currents to establish on long distances. The pro- 
posed mechanism fails therefore to explain the apparent 
poloidal structure of fields around filaments. 



IV. CONCLUSIONS 

We presented a new model for the generation of large 
scale magnetic fields. The mechanism proposed relies 
upon the simple physics of local charge separation which 
occurs naturally in the cosmological context of reionisa- 
tion, at a redshift z ~ 7. As we have seen, the first lu- 
minous objects, such as quasars, emit a radiation which 
ionises the gas filling the Universe. The anisotropy of the 
ionising flux leads to the generation of electrical fields. 
The distribution of the ionising sources together with the 
gas density inhomogeneities imply that the ionising flux 
is inhomogeneous. The fluctuations in the flux generate 
electric currents which in turn induce a magnetic field on 
the scale of the fluctuations. 

With only a couple of simple initial assumptions, 
we could derive interesting properties for the generated 
magnetic fields. In particular, the fields produced in 
our model have amplitudes several orders of magnitude 
higher than in usual plasma physics magnetogenesis mod- 
els. These fields are therefore suitable to serve as seeds 
on large scales for adiabatic collapse amplification up to 
10- 10 Gauss in galaxies. Such large values represent an 
improvement for further amplification by dynamo mech- 
anisms since less time is required in order to account for 
magnetic fields observed in high redshift galaxies. 

Furthermore, relating the magnetic field power spec- 
trum to the matter density power spectrum, we have 
shown that magnetic fields are generated in our model 
preferably on large scales, whereas small scale fields are 
strongly suppressed. This result is true at least in the 
non-saturated regime and we are currently exploring the 
saturated case. Such a property is highly interesting 
again with respect to dynamo amplification theories: the 
small scale fields are weak and therefore their amplifica- 
tion by dynamo represent a less crucial threat. 

We showed also that, in the context of the ISM, our 
mechanism is negligible as a source of magnetic fields 
when compared to usual effects such as dynamo and tur- 
bulence. 

As we mentioned earlier, Subramanian et. al. p3| | and 
Gnedin et. al. |Q explored the generation of magnetic 
field seeds in the context of reionisation. In their model, 
the magnetic field source term is the Battery effect re- 
lying on charge separation. The latter is obtained when 
ionisation fronts either break out from protogalaxies or 
propagate past overdense regions of the Universe. Across 
an ionising front, the temperature varies sharply by sev- 
eral orders of magnitude. The driving physical process 
for charge separation is the thermal pressure of the gas, 
acting in the same way on electrons and ions. Since their 
mass is far smaller than the ion mass, electrons are more 
accelerated by the pressure force than ions and this pro- 
vides charge separation. However, the mechanism we 
presented in this paper relies on a quite different pro- 
cess: The source of charge separation is not differential 
thermal pressure but differential radiation pressure. For 
a particle i = e,p, the acceleration due to thermal pres- 
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sure is a\ oc m i 1 whereas radiation pressure provides an 
acceleration a\ oc m~ 3 . Radiation pressure is therefore 
more efficient by a factor of 



a c/ a p 



3. 10 b 



(42) 



than thermal pressure for accelerating electrons and cre- 
ating electrical fields. This explains why our model leads 
to orders of magnitude higher predictions than in p3[ . 
Thus, the mechanism we described in this paper is likely 
to be dominant in cosmological situations where strong 
anisotropic and inhomogeneous radiation pressure exists. 
Further developments of the model presented in 



this paper, relaxing our geometrical assumptions and 
addressing the magnetic field power spectrum in the 
saturated regime will appear in a forthcoming paper. 
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